Quark charge balance function and hadronization effects in relativistic heavy ion collisions 
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We calculate the charge balance function of the bulk quark system before hadronization and those for the 
directly produced and the final hadron system in high energy heavy ion collisions. We use the covariance 
coefficient to describe the strength of the correlation between the momentum of the quark and that of the anti- 
quark if they are produced in a pair and fix the parameter by comparing the results for hadrons with the available 
data. We study the hadronization effects and decay contributions by comparing the results for hadrons with 
those for the bulk quark system. Our results show that while hadronization via quark combination mechanism 
slightly increases the width of the charge balance functions, it preserves the main features of these functions 
such as the longitudinal boost invariance and scaling properties in rapidity space. The influence from resonance 
decays on the width of the balance function is more significant but it does not destroy its boost invariance and 
scaling properties in rapidity space either. The balance functions in azimuthal direction are also presented. 

PACS numbers: 25.75.Dw, 25.75.Gz, 25.75.Nq, 25.75.-q 

I. INTRODUCTION 

The electric charge balance function for the final state hadrons has been proposed as a probe to study the properties of 
the bulk matter system produced in relativistic heavy ion collisions |ll|-|7|]. Measurements have already carried out both 



nil , we are already able 



in the rapidity space Il8t- ll0ll and in the azimuthal direction llllll . From the data now available 111 Ot 
to see clearly that the charge balance functions for hadrons produced in high energy heavy ion collisions are significantly 
narrower than those for pp collisions at the same energies and they are narrower for central collisions than those for 
peripheral colUsions, indicating a strong local charge compensation in the bulk quark matter system produced in heavy 



ion collisions. The data 11211 further show that the charge balance functions have the longitudinal boost invariance and 
scaling properties in the rapidity space, and these properties hold for either transverse momentum /^^-integrated balance 
functions or those for different pj ranges. 



These features of the experimental data lll0l - ll2ll are rather striking and suggest that such studies should be able to give 



more insights to the understanding of the properties of the bulk quark matter system produced in AA collisions. It is thus 



natural to ask whether such behavior hold also for the quark anti-quark system before hadronization. It is also important 
to see how large the influence from the hadronization and resonance decay. 

In this paper, we propose a simple working model to calculate the charge balance function for the bulk quark anti-quark 
system before hadronization. We introduce the variance coefficient p to describe the local correlation in the momentum 
distribution for the quark and that for the anti-quark if they are produced in a pair The parameter p measures the strength 
of the quark-anti-quark momentum correlation produced in the processes. We study the influence due to hadronization 
process including the contributions due to resonance decay by simulating the hadronization process using a quark combi- 
nation model which describe the final hadron distributions. 

The paper is organized as follows. In Sec. II, we study the charge balance of the quark system before hadronization. In 
Sec. Ill, we study the charge balance function of initial hadron system as well as final hadron system, and compare them 
with that of quark system. Sec. IV gives a brief summary. 

II. CHARGE BALANCE FUNCTION OF THE SYSTEM OF QUARKS AND ANTI-QUARKS 

We recall that , the balance function is in general defined as i^, 

B(A2|Ai) = -[p{b, A.lfl, Ai) - p(fl, A2|fl, Ai) + p(fl, A2IZ7, Ai) - p{b, A2IZ7, Ai)}, (1) 

where p{b,A2\a,Ai) is the conditional probability of observing a particle of type b in bin A2 given the existence of a 

particle of type a in bin Ai . The label a may e.g. refer to all positively charged particles while b refers to all negatively 

charged ones; a may also refer to all particles with strangeness - 1 while b refers to those with + 1 , and so on. For a system 

consisting of many particles, the conditional probability p(i',A2|fl,Ai) is calculated by counting the number A^(/7,A2|fl,Ai) 

of the ab-paiis where a is in bin Ai and b is in bin A2 and the number A^(fl, Ai) of a in bin Ai, i.e., 

,. . , . , N{b,A2\a,Ai) 
p(b, A2\a, Ai) = — .., . , • (2) 

N{a, Ai) 

These numbers can be calculated using the usual two-particle joint momentum distribution function fab{pi,Pi) and single 

particle distribution function faip) or ftip) respectively. They are given by, 

A?(/7,A2|fl,Ai)= { d'pi f d'pifabiPuPi), (3) 

Ja, Jaj 

A^(fl,Ai)= r d'piUpi). (4) 

Jai 

We see that, if a is locally compensated by b, the balance function B(A2|Ai) should have a very narrower distribution. 
In the opposite case, it should be flat. In the case that a is globally compensated by b, e.g., for electric charge balance 



function where a and b denote positively or negatively charged particle respectively, the balance function is normalized to 

unit, i.e. Xa^ B(A2|Ai) = 1. 

A. A working model for the two particle joint momentum distribution functions In the bulk quark matter system 

We consider the bulk quark matter system produced in heavy ion collisions at high energies. We suppose that the 
system is composed of Ng quarks and Ng anti-quarks. We denote the normalized momentum distribution of the quarks 
and anti-quarks by «^(p) and ng{p) respectively. In heavy ion collisions, the bulk matter system consists of new created 
quarks, anti-quarks and the quarks from the incident nuclei. Those quarks from the incident nuclei are referred as the net 
quarks and the new born quarks and anti-quarks are created in pairs. 

To obtain a charge balance function that is narrower than that for the completely uncorrected case, we introduce 
a minimum correlation in the two particle joint momentum distributions in the system. To this end, we construct the 
following working model for the two particle joint momentum distribution for the bulk quark matter system. We assume 
that there is no correlation between the momentum distributions of two different quarks or two anti-quarks. The joint 
distributions are simply the products of the corresponding single particle momentum distributions, i.e. 

UqiiP\^P2) = Ng,Ng,nq^{p{)ng^(p2)(\ - S,n,g,) + N g,(N g, ' 1 )«,,, (/?! K^ (/^l)^^, ,^2 ' C^) 

fq^-qiiPUPl) = N-g^N-g,ng,{px)n-g^_{p2){\ - 6g^,g^) + Ng^iNg, ' l^g^ip l)ng^(p2)S g, ,g, , (6) 

where qi and ^2 denote the flavors of the quarks. For the qq joint momentum distribution, we introduce a correlation 
between the moment distribution of the quark and that of the anti-quark which are produced in the same pair In this case, 
the joint distribution for a quark qi and an anti-quark q2 is given by, 

fqmiPi^Pi) = Ng^Ng,ng^(pi)ng-^(p2) + Ng-, [n'^"g"'{puP2) - ng,(pi)ng-,ip2)]Sg,,g,. (7) 

The single particle momentum distributions are related to «^^"^(pi , P2) by, 

ng(p2) = / d^p 1 n'J'C/'i , Pi), (8) 

nqip) = w.'^qip) + \n„„{p). (9) 

Hence, as long as we know n^-"^{pi,p2) and n„ei{p), we can calculate the two particle joint momentum distributions for 
qq, qq and g'^-system. 

To calculate n'lT {pi,P2), we adopt the picture of the hydrodynamic theory. Here, we assume the local thermalization 



and collectivity in the system 113 



1511 . Hence, in the co-moving frame of the fluid cell, due to local thermalization, we 



take a Boltzmann distribution for the single quark or anti-quark distribution, i.e., 

where the supscript * denote that these quantities are in the co-moving frame, K2 is the Bessel function, m is the mass of 
the constituent quark (340 MeV for u ox d quark and 500 MeV for strange quark), and E* - yjp*^ + nP- is the energy of 



quark; T is the temperature of the system at hadronization (take as T = 165 MeV 1160). 

For the joint momentum distribution of the quark and the anti-quark produced in the same pair, we use the covariance 
coefficient p to describe the correlation between them. We recall that for a joint momentum distribution fof a qq-sysXem., 
the covariance coefficientp is defined as p = cov{pq,pij)/vai{pq), where coy{pq,pq) = {pq ■ p^) - {pq) ■ (pq) and var(pq) = 
(/??) - (pq)^. We take the joint distribution for the g'^-system in the co-moving frame of the ^^-pair in the Cholesky 
factorization form, i.e.. 



pair,, , ,. 1 , ., A'PP'I . , n , PcrPPq , 



The covariance parameter p describe the strength of the correlation. If p = 0, there is no correlation between the momen- 
tum distribution of the quark and that of the anti-quark and we obtain the factorized form. For p very close to unity, we 
get a maximum correlation between the momentum of p^ and pq, where the probability is non-zero only when pq - pq. 
In general -1 < p < 1, and p > means short range compensation of q and q while p < means the opposite. 

The joint distribution n'Tg\Pq,Pzi) in the laboratory frame is obtained from n'^""^*(p*q, pi)- Here, we first make the 
Lorentz transformation (fi) from the co-moving frame of the fluid cell to the laboratory frame to obtain n''"" {pq, Pq,fi), 
then sum up the contributions from different fluid cells in the system with different collective velocities, i.e., 

n'qf'iP,' P^) = J h(fi) ^iTiPq, Pq,P) d'/3, (12) 

where h(fi) is the so-called velocity function which corresponds to the velocity distribution of the fluid cell in the system. 
The velocity function h(fi) is normalized to unity and can be decomposed into the longitudinal part hi and the transverse 
part lij_. The longitudinal velocity yS, is usually replaced by the rapidity y. The azimuthal dependence is isotropic, we inte- 
grate it out and obtain, j liifi) (Pp - J hL(y)h±(f^±) dydf5_L. This velocity function liifi) determines, together with the mo- 
mentum distribution of the quark and anti-quark in the fluid cell, the single quark spectrum thus the inclusive momentum 
distribution of the hadrons after hadronization. In practice, it is parameterized by fitting the data for the hadron momentum 



distributions with the aid of hadronization models. According to the transparency observed in experiments! 17 , and be- 
cause of that the observed rapidity spectra of hadrons show a roughly Gaussian shape in the full rapidity range lllSll . we take 



the transverse part as a uniform distribution between [0,y8™^] and parameterize the longitudinal part in a Gaussian-hke 
form, 

hdy)^^^ -e-M°/-\ (13) 

2cr^r(l + i) 

The free parameters a, <t andyS'""^ are fixed using the data for the rapidity and the pr spectra of hadrons. For example, in 

the following of this paper, we just use the resul ts obtai ned by fitting the data of rapidity and pr spectra of final ha drons i n 

central Au+Au collisions at -Jsnn = 200 GeV 



iisot 



m 



1911 with the aid of the combination model for hadronizationll20 
The results are a = 2.40, a- = 2.54 and/3^"-* = 0.30 for u and d newborn quarks and a = 2.36, cr = 2.73 and/?'™* = 0.34 
for strange quarks. The numbers of light and strange (anti-)quarks and momentum distribution of net-quarks from the 



colliding nuclei have been fixed in Ref . 12211 . 



B. Charge balance function of the bulk quark anti-quark system 

Having the joint momentum distribution functions, we can calculate the charge balance function in a straight forward 
way. In the following, we present the results in rapidity space for diff'erent transverse momentum intervals. In practice, 
the balance function in rapidity space is often rewritten as a function of the rapidity diff'erence 6y = ya - yt between two 
particles in a limited window Jm,, i.e., 

D /x \ \ ^ f ^bai^y' 3'h') - NaaiSy, Jh') Nabi^y, Jh) - Nhbi^y, y„) 1 

J3ab(oy\yw) = T 777 — ^ + ^77 — \ ■ '■^^■' 

Since quarks of different flavors posses difl'erent electric charges, it is not straight forward to extend the definition of the 
the electric charge balance function given by Eq.(l) or (fl4t to the quark anti-quark system. There is no direct extension 
of Eq. (fl4l i to such cases. We have many diff'erent possibilities at the quark level, e.g., 

„(ci),. I , 1 V eaebNah(6y,y^) 

1 y^y-' 2Nf Z^ N,{y,,) 

where both a and b run over all the quarks and the anti-quarks, Nf is the number of flavor involved. We can also defined 

it as, 

„(c3), ^ , , If Sfl.fe sgn(e aeb)Nba(6y, y,,) Zaj, sgn(e aeb)NabiSy, y„)) 

2^ LaNaiyw) LbNbiyw) ' 



where a - u, d or s while b - u, d or s represent the positively and negatively charged particles respectively. We may also 
define the baryon number balance function Bq{Sy\y^i?) for the quark anti-quark system instead, which is given by, 

where the summations over a and b run over all different flavors of quarks and those of anti-quarks, and B„ and Bi, stand 
for the baryon numbers. We can also defined it as, 

„(b2),, , . l(I,a.bBaBbNUSy,yw) lLa.bB„BhNab{Sy,yw)^ ,,_,. 

Bf\5y\y,^ = -- ^ „,^ , ^ . + ^ .,,,,, , > (19) 

2* LaB-„N„(y„) LhBiNbiyu.) ' 

where a denotes all the quarks of different flavors and b all the anti-quarks of different flavors respectively. All these 
definitions satisfy J ddyBqidylyw) - 1 . 

We note that, so far as the kind of correlations between the momentum distributions of the quarks and that of the 
anti-quarks described in the working model presented in Sec. A are concerned, all these definitions do not make much 
differences. More precisely, in the working model presented in Sec. A, only a correlation between the momentum of the 
quark and that of the anti-quark from the same qq pair is introduced as given by Eq. (fTTT i. There is no correlation between 
the quarks and anti-quarks from different pairs and there is no difference between different flavors. In this case, all the 
definitions given by Eqs.( ll5m"9l ) are equivalent in the sense that they are all different suppositions of the correlations given 
by Eq.lfTTT) for different flavors and Eq.lfTTt does not distinguish between different flavors. The only differences come from 
the net quark contributions where no strange quark exists. 

For comparison, we made the calculations using the different definitions Eqs.( ll5m"9] l and the results are indeed similar. 
In the following part of this section, we show the results obtained by using Eq. ( fTTl i. 

We first study the case where p = 0. In this case, there is no correlation between the momentum distribution of the 
quarks and anti-quarks. The balance is obtained only from the global flavor compensation of the new created quarks and 
anti-quarks. This is also the minimum compensation in the produced system. In Fig. [I](a), we show the results of y^^l 
for different rapidity positions with transverse momentum p^ integrated. In Fig. [T](b), we show the results for p = 0.5 
and a comparison of the results for different values of p is given in Fig. \T\{c). 

From the results, we see that in all cases, also for p = 0, the balance function Bg{6y\yK) decrease with increasing 
dy showing a local compensation of the electric charge in the rapidity space for the bulk quark system. It is also clear 
that Bq(dy\yn) decreases faster with increasing 6y for larger value of p indicating stronger local charge compensation. 
We also see that Bqidylyn) does not change much for different rapidity window with the same window size showing the 
longitudinal boost invariance. This is hold for different values of the variance coefficient p. 
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FIG. 1: The electric charge balance function B^(dy\y„) for the bulk quark system for same window size as a function of dy at the 
variance coefficient p = in panel (a) and p = 0.5 in panel (b), respectively; A comparison of the results at different values of p is 
shown in panel (c). 



The existence of the approximate boost invariance for the charge balance function for the bulk quark system can easily 
be understood. We note that by looking at the different rapidity window in the case that the window size is much smaller 
than the total rapidity range of the bulk quark system, we are in fact looking at different fluid cells. Since we do not 
differentiate these fluid cells in any significant way, the results should be similar This results in similar charge balance 
function as indicated by the calculated results shown in Fig. \T\ In other words, the boost invariance of the charge balance 
function just reflects the homogeneity of the fluid cell at hadronization in different rapidity windows. 

We continue to study the dependence of the balance function Bg{6y\yH,) on the window size and/or transverse momen- 
tum. In Fig. |2](a), we show Bq(6y\yH) in the different rapidity positions with the same window size y^: = 1 and in Fig. 
|2](b), we show Bg{6y\yH,) at different window sizes y„ - 1,2,3,4. We see that Bgidylyn) varies with window size and 
becomes flatter with increasing window size. This qualitative feature is naturally expected from the definition since the 
balance function is normalized to unity but the range of the allowed values of 6y becomes larger for the larger window 
size. This effect can be eliminated by scaling the balance function Bg{6y\y„) with the factor 1 - ^y/ly^l as suggested in 
Ref yl], i.e. we study the scaled balance function, 

Bg{Sy\y,^) 



BsiSy) 



(20) 



1 - Sy/y„ 

In Fig. |2jc), we show the results obtained for the scaled Bs{6y) of the bulk quark system. We see clearly that the scaled 
balance functions fall on one curve showing that they are independent of the size and position of rapidity window. For 
comparison, we also present the balance function in the full rapidity region (open cross) for the case that the net charge 
of the system is taken to be zero. We see that the result is also consistent with those for the limited rapidity windows so 
far as the scaled balance function is studied. This is very nice feature since it suggests that the scaled balance function 



8 



for particles in the limited rapidity window can indeed be regarded as an example for the charge balance function of the 
system. 
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FIG. 2: The pr -integrated Bg(Sy\y„) of the constituent quark system at different rapidity positions with same (panel a) and different 
(panel b) window sizes, as well as the B,(6y) (panel c). Correlation coefficient p is taken to be 0.3. 



We emphasize that these properties of the balance functions of the bulk quark system are results of the momentum 
distributions of the quarks and anti-quarks in the system . These distributions including the correlations given by Eq. 
(fTTT i are results of the local thermalization and collectivity for the system produced in relativistic heavy ion collisions in 
the hydrodynamic theory. These qualitative features for the charge balance functions for the bulk quark system before 



hadronization is consistent with those for the final hadrons as observed by STAR Collaboration at RHIC lll2ll . 

In Figl3l we show Bg{Sy\yw) and Bs{Sy) in different rapidity windows and in the different pr ranges. We clearly see 
that the scaling properties of balance function still hold in the different p-p ranges. We can also see that the width of the 
scaled balance function decreases with increasing pr- This is because, in general, the quarks and anti-quarks with larger 
Pt come from the fluid cell with larger transverse flow, which results in a smaller longitudinal rapidity interval and hence 
smaller width for balance function. Such a feature was expected earlier at the hadron level [[3| and observed in central 



Au-hAu collisions at -JJnn = 200 GeV 



III. CHARGE BALANCE FUNCTIONS OF THE HADRON SYSTEM 



With the momentum distribution functions of the bulk quark system discussed in last section, we study the charge 
balance functions of hadrons produced in the hadronization of this system. We compare the results obtained for the 
directly produced hadrons and those of the final state hadrons with those for the quarks and anti-quarks to study the 
influence of the hadronization and resonance decay on the balance functions. 

We describe the hadronization of the bulk quark system with the (re-)combination or coalescence mechanism. Such a 
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FIG. 3: The Bg(6y\y„.) of quarlc system (top panels) at different rapidity positions with different window sizes as well as the S.,((5y) 
(below panels) in the different pr (GeV/c) ranges. Correlation coefficient p is taken to be 0.3. 

hadronization mechanism is tested by various data and is implemented in different forms such as the quark recombination 



model 1123 



2411 . the parton coalescence model 125 



2611 . and the quark combination mode (SDQCM) 120 



211]. All these 



models are tested against the various features of the hadrons produced in heavy ion collisions at high energies. Here, in this 



paper, we use SDQCM ]20, 



2111 for our calculations since this model takes the exclusive description and is implemented 
by a Monte-Carlo program so that can be apply to calculate the balance functions for the directly produced hadrons as 
well as the final hadrons after the resonance decays in a very convenient way. Also, this model guarantees that mesons 
and baryons exhaust all the quarks and anti-quarks in the deconfined color-neutral system at hadronization. 



A. Charge balance functions in rapidity space 



We insert the momentum distributions including the correlations given by Eq. (fTTl i to determine the momenta of 
the quarks and anti-quarks before hadronization. We then apply the quark combination rules as implemented in the 



Monte-Carlo program of SDQCM 112 ill to calculate the momentum distribution of the directly produced hadrons. Those 



resonances will decay accordingly and the momentum distributions are simulated also in the program by using the material 



from the particle data group 1127 1 



In Fig. m we show the results for the p^-integrated balance functions for the directly produced hadrons. Here, in Fig. 
|3a), we see the results in different rapidity windows with the same width y^ - 1, while in Fig. |4tb) and (c), we see the 



10 



results at different window sizes as well as the scaled function Bs{6y). In Fig|5] we show the corresponding results in 
different pr ranges. 
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FIG. 4: The pr-integrated B(6y\y,i) of initial hadron system at different rapidity positions with same (panel a) and different (panel b) 
window sizes, as well as the Bs(6y) (panel c). Correlation coefficient p is taken to be 0.3. 
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FIG. 5: The B(6y\yt,) of initial hadron system (top panels) at different rapidity positions with different window sizes as well as the 
6j((5>') (below panels) in the different pj- (GeV/c) ranges. Correlation coefficient p is taken to be 0.3. 



From these results, we see that both the longitudinal boost invariance and rapidity scaling for the balance functions 
are hold for the hadrons directly produced in the quark combination mechanism, either for the ;97--integrated quantities 
or those for different pr ranges. This is in fact not surprising because the formation of hadrons in this hadronization 
mechanism is realized by the combination of two or three nearest quarks/antiquarks in momentum space. This means that 
the combination happens locally and does not destroy the locality nature of charge balance of the system. 

We further study the resonance decay contributions by calculating the balance functions for the final hadrons where 
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decays of the resonances are taken into account. We show the corresponding results in Figs. |6]and|7] 



^^0.7 

^0.6 

to 

0.4 
0.3 
0.2 
0.1 




(a) 



-0.5<y< 0.5 
0<y<1 
1<y<2 
1.5<y<2.5 



0.5 



0.5<y< 0.5 
1.0<y<1.0 
1.5<y<1.5 
2.0<y< 2.0 



5y 




0.5<y< 0.5 
1.0<y<1.0 
1.5<y<1.5 
2.0<y< 2.0 




FIG. 6: The pr-integrated B(6y\y„-) of final hadron system at different rapidity positions with same (panel a) and different (panel b) 
window sizes, as well as the Bs(6y) (panel c). Correlation coefficient p is taken to be 0.3. 
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FIG. 7: The B{Sy\y„) of final hadron system (top panels) at different rapidity positions with different window sizes as well as the B^iSy) 
(below panels) in the different pr (GeV/c) ranges. Correlation coefficient p is taken to be 0.3. 



From Figs. |6]and|7] we see that both the boost invariance in rapidity space and the scaling property still preserved after 
the contributions from the resonance decays are taken into account. Together with those results given in Figs. |4]and|5] 
these results show clearly, although there are definitely influences from hadronization and resonance decay on the form 
of the charge balance functions, these effects do not significantly influence the boost invariance and the scaling in rapidity 
space. 

The influences from hadronization and resonance decay to the balance function can be studied more quantitatively by 
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FIG. 8: In (a), we see the averaged widths {Srf) of the balance functions for the quark and anti-quark, the directly produced and the final 
hadron systems as functions of p. In (b), we see the averaged widths {Srj) of the balance functions for the daughter hadrons from all 
the resonance decays, those for the daughter hadrons from the meson decays and those from the baryon decays, respectively. The band 
area represents the experimental of (^77) for charged particles in central Au+Au ^snn = 200 GeV 11211 . 

calculating the averaged width of the balance function, which is defined as, 

jj'" B{6y\y„) 5y dSy 



(Sy)^ 



£'" B{5y\y,)d5y 



(21) 



We note that the averaged width {5y) is in general a charactering quantity describing the radius of charge balance of 
the system. For the final hadron system in heavy ion collisions, it can be sensitive to different effects such as delayed 
hadronization or hadron freeze-out ill 



021 



3011 . multiplicity effect 



m. 



2811 . possibly highly localized charge balance at freeze-out 112911 . transverse flow 
32II and hadronic weak decay. Here, by comparing the results for the balance functions of 
the quark anti-quark system with those for the initial hadrons and those for the final hadrons, we can study the magnitudes 
of the influences of hadronization and those from resonance decay. 

In Fig. [8ja), we show the results for the averaged width (^77) of the balance function for the directly produced hadrons 
as the function of p, compared with that for the quarks and anti-quarks and that for the final hadrons. Here, we show 



the results for pseudo-rapidity rj in order to compare withjhe experimental data now available lllll . We choose also the 



same pseudo-rapidity and pj regions as the experiments llllll . i.e. I77I < 1, 0.1 < 



In Fig. [8ja), the data of {Stj) in central Au+Au collisions at yfstm - 200 GeV 
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< 2.0 and 0.2 < pr < 2.0 GeV/c. 
] is shown as a band area. We see 



clearly that, in all the three cases, the averaged width {drf) decreases with the increasing p. We see also that, the {drf) for 
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the directly produced hadron system decreases with the increasing p in exactly the same way as that for the bulk quark 
system does. The difference between {Srj) for the directly produced hadron system and that for the bulk quark system is 
almost a constant 0.04 for all the different values of p. This is because, as mentioned above, the combination of quark(s) 
and/or anti-quark(s) in neighbor does not change the electric charge balance in any essential way. However, the formation 
of electrically neutral hadrons may delay the charge balance in momentum space. We take a quark and an anti-quark 
from a given q'^-pair as an example. As given by Eq.lfTTTi. their momentum distributions posses a correlation measured 
by the covariance coefficient p. If both of them enter into the respectively charged hadrons in the combination process, 
the correlation will pass to the hadronic level. However, if one of them enters into an electrically neutral hadron, the 
correlation will be lost in the charged hadrons. This will decrease the local charge balance at the hadron level. 

From Fig. [8ja), we also see that the resonance decay contributions change {5t]) significantly. It is also interesting to 
see that these contributions strengthen the local charge balance for relatively small values of p but weaken the balance 
for larger values of p. This indicates that the decay contributions dilute the balance functions quite significantly. To see 
where these different behaviors come from, we calculate the averaged width {6y) for those hadrons from resonance decay 
separately. We note that the influence of resonance decay to the charge balance function is in general different for hyperon 
decay from that for vector meson decay. The decay of the hyperons such as A ^ p/r and S" -^ Kn produces a pair of 
charged daughter particles with quite narrow rapidity interval, e.g. about one third for A — > p/r and E*' — > Att, due to the 
small amount of energy released in the decay process. This leads also to a smaller {5t]) for the charge balance function. 
However, in vector meson decay such as p° — > n^n' and K*^ — » K^n, the energy released is much larger leading to a 
much larger rapidity difference between the daughter particles, e.g. up to 1.7 forp" -^ n^n^ and 1.3 for A'**' -^ K^n' in 
the rest frame of parent particle. To study this effect in a more quantitative manner, we calculate the averaged width {drj) 
of the balance function only for the daughter particles from baryon or meson decay respectively. The results are shown in 
Fig. [8](b). Here, we see clearly that, the averaged width {dr/} for the daughter particles from baryon decay is indeed much 
smaller than those from meson decay. We see also that the charge balance for the daughter particles from baryon decays 
is dominated by the decay effect which leads to an averaged width of about one third. However, for those from meson 
decays, the charge balance is dominated by the effect from the mother particles. 

B. Charge balance In the azimuthal direction 

The charge balance in the azimuthal direction for hadrons in high energy heavy ion collisions can be sensitive to jet 
production. Experimental studies have already been carried out by STAR Collaboration for hadrons of different pr regions 
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FIG. 9: Balance function Ba-i(6if>) of the quark anti-quark system (a), tlie directly produced hadron system (b) and the final hadron 
system (c) as functions of dip in the pseudo-rapidity region - 1 < tj < 1 and 0.2 < pj < 2.0GeV/c. A comparison of the results for 
the final hadron system at p = 0.5, 0.6 and 0.7 and the experimental data for all charged particles with 0.2 < pj < 2.0GeV/c in central 

■.Im is ] 



Au+Au collisions at -Js^n = 200 GeV in Ref.( 



; given in (d). 



I nil . It is thus also interesting to see how the charge balance function behaves for the bulk quark matter system and the 
resulting hadrons. 

The balance function of hadrons in the azimuthal direction is defined similarly to that in rapidity, 



D ,^^ ^, I (Nhai6<f>,lP)-Naa{6lP,(f>) NalAS<P, Cp) - Nhb{64>, <f>) 

B,.,.-,m,^) = - I j^^ . j^^ 



(22) 



where, e.g. the quantity Nha(6(/), (p) is defined as the number of pairs of particles with the particle a flying at an angle (p 
(measured with respect to the reaction plane) and the particle b at an angle between and + 6(f). In the following we 
study the azimuthally averaged balance function 



Bab,a:l(S4>) 



pin 



d4>Bat(64>,(f>). 



(23) 



Having the Monte-Carlo program at hand, the extension of the calculations mentioned above to azimuthal direction is 
straight forward. We show the results obtained for the quark, the directly produced hadron and the final hadron system at 



different p values in Fig. |9] Comparison with the available data llllll is also given in the figure. 

From the results for the bulk quark system Fig. |9ja), we see clearly that the dependence of the quark charge balance 
function on the variance parameter p is quite obvious. For p close to unity, the momentum of the quark and that of the 
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anti-quark produced in pair are closely correlated, and we see a sharp peak at 64> - 0. For p - Q, there is no correlation 
and the balance function is almost a flat function showing only the influence from the global charge compensation. 

Comparing the results in Fig. |9jb) with those in Fig. |9ja), we see that the charge balance functions in the azimuthal 
direction for the directly produced hadron system are slightly broader than the corresponding results for the quark system, 
showing a slightly loose correlation. This is similar to the case in rapidity space studied in last subsection. However, the 
influences from the resonance decays are quite significant in the azimuthal direction. We see quite significant differences 
between the results for the final hadrons and the corresponding results for the hadron system before resonance decay. We 
see in particular that the very much pronounced peak at 5</i = is smoothed by the decay influences. This is also obvious 
since such strong correlation can be destroyed by the resonance decay because of the transverse momentum conservation 



in the decay processes. From Fig. |9jc), we see that the data llllll is well be described except for the peak at 6(p — 0. This 



peak could be an indication of jet contribution which is not included in our calculations. 

IV. SUMMARY 

In summary, we have calculated the charge balance functions of the bulk quark system before hadronization, those for 
the directly produced and the final hadron system in relativistic heavy ion collisions. The momentum distributions for 
the quarks and the anti-quarks in the bulk quark system are taken as determined in the hydrodynamic picture with local 
thermalization and collectivity. A correlation between the momentum distribution of the quark and that of the anti-quark 
is introduced if they are from the same new produced qq pair and the correlation strength is described by the variance 
coeflicient p. Our results show that the charge balance functions for the bulk quark system have the longitudinal boost 
invariance and the scaling behavior in rapidity space. Such properties are preserved by the subsequent hadronization via 
combination mechanism and the resonance decay, although both hadronization and resonance decay can influence the 
width of the balance function. With the same inputs, we also studied the balance function in the azimuthal direction. 
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